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Q ■ ABSTRACT 

%^ I This work aims to study the distribution of luminous and dark matter in Coma early- 

^ . type galaxies. Dynamical masses obtained under the assumption that mass follows 

' light do not match with the masses of strong gravitational lens systems of similar 

velocity dispersions. Instead, dynamical fits with dark matter halos are in good agree- 
ment with lensing results. We derive mass-to-light ratios of the stellar populations 
from Lick absorption line indices, reproducing well the observed galaxy colours. Even 
, in dynamical models with dark matter halos the amount of mass that follows the light 

increases more rapidly with galaxy velocity dispersion than expected for a constant 
■r^lj- I stellar initial mass function (IMF). While galaxies around (Toff ~ 200km/s are consis- 

ts . tent with a Kroupa IMF, the same IMF underpredicts luminous dynamical masses of 

galaxies with ctcS ~ 300km/s by a factor of two and more. A systematic variation of 
I the stellar IMF with galaxy velocity dispersion could explain this trend with a Salpeter 

IMF for the most massive galaxies. If the IMF is instead constant, then some of the 
dark matter in high velocity dispersion galaxies must follow a spatial distribution very 
similar to that of the light. A combination of both, a varying IMF and a component 
of dark matter that follows the light is possible as well. For a subsample of galaxies 
I with old stellar populations we show that the tilt in the fundamental plane can be ex- 

^ . plained by systematic variations of the total (stellar + dark) mass inside the effective 

^ • radius. We tested commonly used mass estimator formulae, finding them accurate at 

the 20 - 30 % level. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: kinematics and dynamics 
— galaxies: structure 
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1 INTRODUCTION 

The masses of galaxies are revealed by the gravitational 
interaction of their matter constituents, e.g. by stellar or 
gas kinematics or gravitational lensing effects. The caveat is 
that these effects rigorously constrain only the total amount 
of gravitating mass. The decomposition into luminous and 
dark matter relies on further assumptions. For example, in 
dynamical studies of early-type galaxies it is commonly as- 
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sumed that the stellar mass distribution follows the light. 
Any radial increase in the mass-to-light ratio - if needed to 
explain the observed velocities of the stars - is attributed to 
dark matter. 

The limitation of this approach becomes clear if one 
imagines a galaxy dark matter halo that follows the light 
distribution exactly. Discriminating between luminous and 
dark matter according to their potentially different radial 
distributions would fail in this case. In fact, there would be 
no direct way to dynamically unravel the relative contribu- 
tions of luminous and dark matter to the galaxy mass. It 
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is unlikely that a real galaxy halo follows the light distri- 
bution exactly, yet it exemplifies the intrinsic degeneracies 
in any mass decomposition. On top of this, even the most 
advanced present-day dynamical modelling techniques rely 
on symmetry assumptions and, often, also on the assump- 
tion of a steady-state dynamical system. If the symmetry 
assumptions are strongly violated, dynamica l mass-to-light 
ratios can be biased by a factor of up to two (|Thomas et al.l 
l2007al ). 

To crosscheck the validity of the assumptions in dy- 
namical modelling it is important to compare the resulting 
masses with other independent methods. Total masses can 
be most directly compared to results from gravitational lens- 
ing, which is the first goal of this paper. 

An examination of the mass decomposition requires the 
investigation of a galaxy's stellar population. Since the latter 
provides an independent measure of the stellar mass-to-light 
ratio, the comparison of dynamical and stellar-population 
masses yields further constraints on the dark matter content. 
In the above mentioned case, for example, the dynamical 
mass-to- light ratio T,.dyn would exceed the corresponding 
stellar population value, indicating that some fraction of the 
galaxy's mass is actually dark matter. 

Unfortunately, stellar population models do not pro- 
vide unique stellar mass-to-light ratios. They suffer from 
an incomplete knowledge of the initial stellar mass function 
(IMF), as well as age-metallicity degeneracies. Observations 
in and around our own Galaxy indi c ate that the I MF slope 
flattens below 0.5 M© l|Scald Il986l : iKroupal [iooH '). Recent 
spectroscopic observations of massive ellipticals in the near- 
infrared point towards the low-mass slope of the IMF i n 
these galaxies being steeper (|van Dokkum fc Conrovll2010l ). 
Yet, until now the IMF in distant galaxies with unresolved 
stellar populations remains largely uncertain. This trans- 
lates into a signiflcant indeterminacy of population mass-to- 
light ratios: the steeper the slope at the low-mass end, the 
higher the population mass-to-light ratio. Hence, the stellar 
population approach is not directly conclusive as a probe of 
the mass decomposition in dynamical models. 

Conversely, if dynamical stellar mass determinations 
were free of ambiguities with respect to a dark matter con- 
tamination then they could serve as a measure for the slope 
of the IMF in distant galaxies. The method would be to 
tweak the IMF in the stellar population models until agree- 
ment with the dynamical stellar masses is achieved. 

The ambiguities in both, dynamical stellar masses as 
well as stellar population models, make neither of the ap- 
proaches directly applicable. Nevertheless, comparing dy- 
namical with stellar population models is important to (1) 
narrow down the dynamically plausible range of possible 
IMFs and to (2) delimit the range of dark matter fractions 
compatible with observed stellar population properties. This 
is the second goal of the present paper in which we compare 
dynamical and stellar population masses in a sample of 16 
Coma early-type galaxies. 

Our dynamical models account for both the full variety 
of possible orbit conflgurations in axisymmetric, flattened 
galaxies as well a s for dark matter . In thi s respect we extend 
previous studies. rCappellari et al.l l|2006l 'l. for example, used 
a similar modelling technique to compare dynamical and 
stellar population masses in the SAURON sample, but they 
did not consider dark matter explicitly in their dynamical 



models. The justification for this neglect was the expected 
insignificance of dark matter in the central galaxy regions 
probed by the SAURON observations (robs, max < ^cff ). How- 
ever, measuring only the sum of luminous and dark mass 
makes the compar ison with stellar popula tion models po- 
tentially uncertain. iNapolitano et al.l ()201[)^ analysed a large 
sample of early-type galaxies taking into account the contri- 
bution from dark matter, but their models do not account 
for galaxy fiattening and rotation. 

A subsamplc of the Coma galaxies was recently anal- 
ysed bylCrillo & Gobat (2010). While they used multi-band 
photometry to derive stellar population properties, our ap- 
proach is to measure ages, metallicities and [a/Fe] ratios 
from Lick indices to reduce potential biases in population 
parameters . Like wise, the analysis of SLAGS galaxies by 
iTreu et all (|2010D . combining constraints from gravitational 
lensing and stellar dynamics, relied on multi-band photom- 
etry for the stellar population part. 

Mass-to-light ratios of early-type galaxies are of par- 
ticular interest to understand the tilt of the fundamen- 
tal plane. Virial relations imply that the effective surface 
brightness {I)cS, the effective radius r^s and the central 
velocity dispersion ao in hot stellar systems are not inde- 
pendent of each other. This is revealed by the fundamen- 
tal plane of early-ty pe galaxies (Diorgovski & Davis 1983; 
iDressler et al.|[l987l ). Yet, the observed fundamental plane 
is tilted with respect to the simple case of a virialised, ho- 
mologous family of dynamical systems. This tilt can reflect 
(1) systemat ic variation s of the lumino sity distribution (e.g. 
ISaglia et aL|[l993 : Truiillo et al.ll2004). (2) svstemat ic vari- 
ations of the orbital structure (e.g. lCiotti et al.]|l996l ) or (3) 
systematic variations of the mass-to-light ratio, as a result of 
vary ing stellar populations and/or dark matter distributions 
(e.g. lRenzini fc Ciotti|[l993l ) . Understanding these variations 
allows a deep er insight into the f ormation process of early- 
type galaxies (|Bender et al.lll993 ). 

Most of the above mentioned effects can be factored 
out if additional information about the stellar populations 
and/or the mass distributions are available. Aperture spec- 
troscopy is one way to measure stellar population proper- 
ties. By assuming simple scaling laws it can also provide 
estimates for dynamical masses, such that the contributions 
of stellar population and dynamical effe cts on the funda- 
mental plane tilt can be dis entangled (e.g. lHvde fc Bernardil 
|2009| :Ig raves fc Fabei|[2O10l ). More reliable constraints come 
from radially resolved spectroscopy and detailed dynamica l 
(or lensing) mode ls of galaxies (e.g. ICappellari et al. I I2OO6I : 
iBolton et al.|[2007l ). The third goal of this paper is to follow 
the latter approach and to use the specific information con- 
tained in our two-component dynamical models for further 
investigations upon the origin of the fundamental plane tilt. 

The paper is organised as follows. Sec. [2] reviews the 
galaxy sample and models, in Sec. [3] we compare projected 
masses from dynamical models and from gravitational lens- 
ing against each other. Sec. |4] deals with the comparison 
of luminous dynamical and stellar population masses. The 
dark matter distribution is analysed in Sec. [S] and implica- 
tions for the tilt of the fundamental plane are addressed in 
Sec. [6] The paper is summarised in Sec. [T] 
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Figure 1. Stellar population mass-to-light ratios Tsaip (Salpeter IMF) as function of radius along the major-axis (filled circles) 
and minor-axis (open circles). Note that for some galaxies there are systematic differences between the two sides of a slit which 
are slightly larger than the statistical errors (e.g. along the major-axis of GMP0144). The dotted line in each panel is the light- 
weighted average of Tg^ip within r < r^ff, the solid line corresponds to the stellar mass-to-light ratio T, (jy^ from dynamical 
models. From top-left to bottom-right galaxies are plotted in order of decreasing T, jy^. 



2 GALAXY SAMPLE AND MODELLING 

The sample analysed in this paper comprises 16 Coma early- 
type galaxies in the luminosity range between Mb ~ —19.88 
and Mb = —22.26 (eight giant ellipticals, eight lenticu- 
lar/intermediate type g alaxies). It is almost identical to 
the sample presented in I Thomas et al.l (|2009bl ) . Only the 
galaxy GMP3958 has been omitted for its strong gas emis- 
sion, which hampers reliable stellar population modelling. 

For the dynamical analysis of each galaxy a composite 
of ground-based and HST photometry has been used. Stel- 
lar absorption line data for the kinematics and Lick indices 
come from various long-slit spectra, at least along the major 
and minor axes, but in many cases covering other position 
angles as well. The spectra extend to 1 — 4reff. Details about 
th e photometric and sp e ctroscopic data have b een published 



iMehlert et all (l200(]|).IWegner et al.l (|2002l ). ICorsini et all 
l|2008l ). and lThomas et al.l (|2009bl ). 



2.1 Dynamical modelling 

To the kinematic and photometric data we applied 
our implementation of Schwarzschild's orbit superposition 



techni que llSchwarzschild 1979|) fo r axisymmetric poten- 
tials ([Richstone fc Tremaind 1 198^ : ICebhardt et al. I l2003l : 
iThomas et al.l 12004 l2005al). A d e tailed description of the 
models is given in iThomas et ai] (|2007bl ). Most important 
for this paper are the assumptions about the mass structure. 
We will consider two sets of models. For the first set it is 
assumed that all the mass follows the light, i.e. 



Ty.v 



(1) 



where v is the three-dimensional luminosity density. By con- 
struction, the mass-to-light ratio T here includes the contri- 
bution of both the stellar and dark mass of a galaxy. It is not 
known in advance and obtained by a ^'^-minimisation with 
respect to the kinematical observations. The best-fit T of 
one-component (i.e. self consistent) models will be referred 
to as T,,sc in the remainder of this paper and can be found 
in Tab. [T] 

For the second, and standard set of models we assume 
a mass density of the form 



p = T X + pdm. 



(2) 



The first component again follows the light, while the sec- 
ond one, pdm, accounts for dark matter. Eq.[2]is designed to 
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separate the contributions of luminous and dark matter to 
the total mass of a galaxy. Our basic assumption is that the 
best-fit T of two-component models represents only the stel- 
lar mass of a galaxy and we will refer to it as the dynamical 
stellar mass-to-light ratio T»,dyn in the following. Strictly 
speaking, T.^dyn measures all the mass which follows the 
light distribution, be it stellar or be it dark matter. In this 
respect T,,dyn provides only an upper limit for the stellar 
mass. As soon as there are other mass components which - 
for whatever reason - follow the light, they will contribute 
to T,,dyn as well and the actual galaxy stellar mass-to-light 
ratio will be smaller than T.^dyn- A more detailed discussion 
upon this issue will be given in Sees. 14.41 and 15.21 

According to equation ((2)1, the cumulative (spherical) 
dark matter fraction /oM.dyn of the models inside radius r 
is 



DM.dyn 



(r) 



DM,dyn 



(r) 



MDM,dyn(r') 



(3) 



dyn(r) T*,dyn L{r) -\- MmA,dyn{r) 

where Mtot.dyn, AfoM.dyn and L stand for the cumulative 
total mass, dark mass and light inside radius r. Tab. [T] lists 
the observationally derived dark matter fractions. Ifere we 
consider mostly the region inside rcff , where the average dark 
matter fraction is (/oM.dyn) = 23± 17 %, but our data reach 
out to radii w here the stellar and t he dark matter density 
become equal (|Thomas et al.ll2007bl) . 

We probed two dark matter descriptions. Firstly, log- 
arithmic halos with a constant-density co re of size rn and 
secondly NFW halos (|Navarro et al.lll996l ). The latter pro- 
vide good fits to cosmological A'^-body simulations and have 
a central logarithmic slope of —1. Neither of the two profiles 
includes baryonic halo contraction explicitly (cf. the discus- 
sion in Sec. I5.2|l . Coma galaxies are in most cases better fit 
with logarithmic halos, but the significance over NFW halo 
profiles is marginal. The steeper central slope of NFW ha- 
los implies a higher central dark matter density. But even 
in cases where NFW halos fit better, the innermost galaxy 
regions are still dominated b y the mass-component that fol- 
lows the light (for details see lThomas et al. I l2007bl ). The as- 
sumptions about the halo-density profile have therefore lit- 
tle infiuence on the best-fit stellar mass-to-light ratio T,,dyn 
(typically, T^^dyn from logarithmic or NFW halos differ by 
no more than ATt^dyn ~ 0.5; see also Sec. I4.2|l . 

The gravitational potential of the galaxies is assumed 
to be axisymmetric. Contrasting other methods like using 
Jeans equations, the Schwarzschild technique allows the ex- 
ploration of all possible orbit configurations. The Coma 
sample is unique in being the only larger sample of ax- 
isymmetric dynamical models including dark matter. Previ- 
ous modelling attem pts either assumed spherical symmetry 
llGerhard et al.ll200ll) or did not include dark matter explic- 
itlv (ICappellari et al.|[200^ ). 



2.2 Stellar population models 

Stellar ages, metallicities, [a/Fe] ratios and _R-band stel- 
lar mass-to-light ratios are de termined by fitting th e sin- 
gle stellar population models of Maraston (Il998l. 20051 ) with 
Q-elements overabundance of iThomas et al.l 1 20031 ) to the 
Lick indices H/3, (Fe), [MgFe] and Mgfo. Two initial-stellar- 
mass functions are considered. Firstly, the Salpeter IMF 
(Tsaip, with mass limits of 0.1 M0 and 100 M0) and, sec- 
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observed SDSS colours 

Figure 2. Observed versus predicted colours for the subsample 
of Coma galaxies with SDSS photometry. Colours as indicated in 
each panel. Model predictions are plotted along the vertical axis, 
SDSS observations horizontally. All axes are in magnitudes. The 
mean colour difference A and rms-scatter between models and 
observations (in magnitudes) is quoted in each panel. 



ondly, the Kroupa IMF (Tkiou, with the same mass limits, 
but a shallower slope for stars below 0.5 Mq). The Salpeter 
IMF implies more low-mass stars and a higher mass-to-light 
ratio. In the _R-ban d the scaling between the t wo cases is 
Tsaip ~ 1.56TKrou (|Salpeted[l955l : lKroupall200lf ) . When the 
IMF does not need to be specified, we will refer to stellar- 
population mass-to- light ratios as T^sp- 

Stellar population properties are calculated at each ra- 
dius with observations. Salpeter 7?-band mass-to-light ra- 
tios are shown in Fig. [1] The light-weighted averages within 
rcff are indicated by the dashed lines and form the basis 
for the remainder of this paper (Tab.[T]). For some galaxies 
there are systematic differences between the two sides of a 
slit which are slightly larger than the statistical errors (e.g. 
along the major-axis of GMP0144). These systematic un- 
certainties are not included in the errors quoted in Tab. [1] 
Irrespect ive of the metalicity gradients present in early-type 
galaxies (|Mehlert et al.|[2003l ). mass-to-light ratio gradients 
in the R-band are generally small for the Coma galaxies (cf. 
Fig. [T]). The results presented here do not depend signifi- 
cantly on the averaging radius. Its choice is driven by the 
most massive sample objects. Their spectroscopic data reach 
only out to robs, max ~ reff and, for the purpose of homogene- 
ity, we restrict the averaging in other galaxies to the same 
radius, even if the data extend further out. 

For a subsample of the Coma galaxies studied here, 
multi-b and photometry from the Sloan Digital Sky Survey 
CSDSS: [York et al.|[200ol ) is available. Fig. [2] compares ob- 
served SDSS colours with the predictions of our best-fit SSP 
models. We have applied the same colour corrections to the 
models as discussed in .Saglia et al. (201 0) plus an additional 
i — z = —0.05 (Maraston, private communication), to take 
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galaxy "'"♦.dyn Tf*,sc /DM,dyn C^eff Tkfou "^Salp 

GMP NGC/IC [Uq/Lr^q] [M0/Lfl,0] [km/s] [Mq/L^.q] [M0/Lfl,0] 



(1) 


(2) 


(3) 


(4) 




(5) 




(6) 






(7) 




(8) 


m /I /I 




A c;n -i- A 

4.0U ± U.DU 


V AA 
I .UU 


A 
U 


dn+o.io 

■^^-0.05 




o it U 




Q 
O 


AQ + A 1 V 
Uo It U. i / 


A 


V/l + A Ofi 
/ 4 It U.^D 


noQO 




nn + A c;a 
o.UU It U.oU 


a c;a 
D.OU 


A 
U 


ncr + O.lS 
^"^-0.05 




z it u 





Q 
O 


+ A OO 
O / It Ki.ZZ 


r; 
O 


Oo It U.o4 


UiOu 




O ^^A + A OA 
z.DU ± U.ZU 


Q AA 
O.VU 


A 
U 


-,1+0.05 
^^-0.05 


1 oq 


Z it U 


Q 
O 


1 
1 


/I A + A Afi 
4U It U.UD 


Z 


1 O + A 1 A 
iy It U. iU 


1 1 '7« 




O OA J- A OA 

Z.UU ± u.zu 


O A 


A 
U 


,1 +0.05 
■^^-0.04 


OAC^ 


o it U 


Q 
O 


-1 
1 


-1- A A/I 

04 it U.U4 




/1 1 -1- A Afi 
4i it U.UD 


1750 


4926 


6.00 ±0.75 


7.00 





-ic:+0.47 
^''—0 05 


279 


1 ± 1 


6 


3 


41 ±0.50 


5 


33 ±0.79 


1990 


IC 843 


10.00 ± 1.00 


10.00 





Q^ + 0.01 

"^-0.01 


281 


3±1 


1 


3 


03 ± 0.30 


4 


74 ± 0.47 


2417 


4908 


8.00 ±0.75 


8.50 





nK+0.20 

'^^-0.04 


211 


0±1 


7 


3 


78 ± 0.43 


5 


91 ±0.67 


2440 


IC 4045 


6.50 ±0.50 


7.00 





A7+0.02 
-0.02 


225 


9±0 


9 


3 


46 ± 0.27 


5 


41 ±0.42 


3414 


4871 


4.00 ± 0.62 


6.00 





46+0-06 
^"-0.27 


169 


6±1 


3 


3 


20 ±0.29 


4 


99 ±0.45 


3510 


4869 


5.50 ± 0.50 


6.00 





"^-0.05 


177 


7± 1 


7 


3 


49 ± 0.65 


5 


45 ±1.02 


3792 


4860 


8.00 ± 1.00 


9.00 





iq+0.24 
^•^-0.05 


284 


0±1 


9 


4 


27 ±0.31 


6 


66 ±0.48 


4822 


4841A 


5.50 ± 1.00 


6.50 





CI +0.14 
^^-0.27 


272 


1 ±2 


7 


3 


24 ±0.33 


5 


07 ±0.52 


4928 


4839 


8.50 ±2.00 


10.00 





no+0.35 
■^^-0.21 


314 


8±2 


9 


3 


34 ±0.46 


5 


22 ±0.72 


5279 


4827 


6.50 ± 0.50 


7.00 





1Q+0.21 
^"-0.07 


244 


1 ± 1 


2 


3 


16 ±0.45 


4 


93 ±0.71 


5568 


4816 


6.00 ± 1.00 


7.00 





cq+0.15 
•^•^-0.25 


233 


4±1 


7 


3 


91 ±0.33 


6 


11 ±0.52 


5975 


4807 


3.00 ±0.50 


4.00 





9Q+0.05 

^■'-o.oi 


195 


9±0 


8 


3 


07 ±0.24 


4 


80 ±0.37 



Table 1. Dynamical parameters of Coma galaxies. (1, 2) galaxy identification (GMP from [Godwin et al.|[l983l) . (3) best-fit dy- 
namical dyn (R-band) in models that explicitly account for dark matter (the quoted errors include all models that devia te by 
less th an Ax'^ 1 from the best-fit model. For a more detailed discussion of the errors the reader is referred to Thomas_e^_alJ 
2005al .'l (4) best-fit dynamical T,,sc (-R-band) assuming that all the mass follows the light. (According to Fig. 2 in lThomas et al.l 
2007bl the formal errors on Tt^sc are smaller than those of jy^. However, since the assumption that mass-foUows-light does 
not yield satisfactory fits to the kinematics we don't give errors for T«^ac here.) (5) dark matter fraction /DM,dyn within r^ff. (6) 
galaxy velocity dispersion a^ff (inside r^g). (7, 8) stellar population mass-to-light ratios for the Kroupa IMF (Tkiou! R-band) 
and Salpeter IMF (Tg^ip; i?-band), respectively. The mass- to-light ratios are light- weighted averages within r^fj . Only radii with a 
stellar-population age r ^ 14 Gyr are considered. 



into account the recent improve ments in the cal ibration of 
the Maraston SSP models (Mar aston et al.ll2009l ). The ver- 
tical error-bars indicate the 68% confidence region derived 
from the observational errors. Model colours are averaged 
inside r^f!. They fit well to the SDSS colours with average 
diff erences smaller t h an < 0.01 mag. 

iGrillo fc Go"ba3 (|2010l ) used observed SDSS colours to 
derive photometric stellar population parameters for some 
of our Coma galaxies. In contrast to the analysis presented 
here, they (1) assume a solar metallicity for all galaxies (but 
see the middle panel of Fig. Illf) (2) allow for an ext ended 
star-formation history and (3) use the lMarastonI (|2005l ) mod- 
els without colour corrections. In Fig. [3] we plot their pho- 
tometric T^f'o^ (scaled to the Kroupa IMF) against our 
Txrou- On average, both approaches yield consistent results 
{{'^KTou/TfSou) = 1-11), though the rms-scatter (±0.35) is 
large. 



3 PROJECTED MASSES FROM DYNAMICS 
AND LENSING 

As it has been stated in Sec. [T] the dynamical models rely 
on the assumption of axial symmetry and steady state dy- 
namics. To check how accurately these assumptions are ful- 
filled in real galaxies we first compare our dynamical mod- 
els against gravitational lensing results. The latter constrain 
the total projected mass inside a cylinder delimited by the 
Einstein radius rEin of the lens and a re less affected by sym- 
metry assumptions (|Kochaneklli99ll ). 



3.1 Lens selection 

The Einstein radius of a gravitational lens results from two 
independent properties of a lensing configuration. Firstly, 
from the physical deflection angle that the lensing galaxy 
gives rise to according to its gravity. It depends only on the 
mass distribution of the foreground galaxy. Secondly, from 
projection factors that depend on the distances of the fore- 
ground lens and the background source, respectively. For the 
Coma galaxies we only know their mass distributions, but 
they are not part of real lenses. To compare them with ob- 
served gravitational lenses, we need to define an appropriate 
lensing configuration for each Coma galaxy. Here we do this 
implicitly by seeking for lensing galaxies that accidentally 
happen to fall on a linear relation rEin(o"cff) between the 
Einstein radius and the effective velocity dispersion. Defin- 
ing a fiducial Einstein radius for each Coma galaxy accord- 
ing to the same rEin(o"cfi) then ensures, that there is at least 
a subsample of real lenses with similar lens configurations 
at a given (Tcfi. The Coma galaxies can be compared to the 
lenses in such a subsample, but not to the rest of the lensing 
galaxies. The form of the selection function is arbitrary, any 
other function would serve equally well. We choose a linear 
relation for simplicity. 

Fig. [4] s hows th e observed Einstein radii rEin of SLACS 
lenses from lAuger et al. ( 2009) against their velocity dis- 
persion^. The large circles show two different subsamples 



Note that for Fig. |4] we have used the correction of 
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Figure 3. Stellar population Yxrou (K roupa IMF) from SDSS 
colours (x-axis; cf. iGrillo fc Gobatllioiol) and from Lick indices 
(y-axis; cf. Tab.[TJ- 

of lenses constructed as outlined above. The lenses are se- 
lected to deviate by less than 0.2 kpc from two arbitrary 
linear selection functions rEin(croff) = 0.025 (Tcft — 3.0 and 
'"Bin(o"efi) = 0.025 (Tcff — 1.9. After having selected the lenses 
we fit a straight line to each subsample in order to determine 
the actual best-fit rEin(o"cff) that is used to define fiducial 
Einstein radii for Coma galaxies. The fitted relations dif- 
fer only slightly from the original selection functions. For 
subsample SL0.5 (dotted) we get 

^ = 0.0238 X - 2.6709 (4) 

kpc km/s 

and for subsample SL0.75 (solid) 

^ = 0.0249 X - 1.8564. (5) 

kpc km/s 

The two subsamples are constructed as a compromise be- 
tween (1) ending up with a sizeable number of galaxies in 
each subsample and (2) yielding sufficiently different sub- 
samples to allow for a comparison between Coma and lens- 
ing galaxies at different physical scales. For subsample SL0.5 
we get an average (r-Ein/?'cft) ~ 0.5 (11 lenses), while for sub- 
sample SL0.75 it is (rEin/rcff ) ~ 0.75 (17 lenses). 

3.2 The total mass 

The large symbols in Fig. [5] show projected, integrated 
masses of Coma galaxies and SLACS lenses as a function 
of galaxy velocity dispersion (JcEf. Coma masses are calcu- 
lated from the integral 

/ 10 '-off /-air (•'■EinCo'ctt) 

Az j dip prdr, (6) 

-10r„ff Jo Jo 

ICappellari et al to transform the measured SDSS aper- 

ture dispersions into cTgff. 



Figure 4. Einstein radii r Ein and effective ve locity dispersions 
(Teff of SLACS galaxies from lAuger et al] |200^ are shown by the 
crosses. Dotted open circles highlight lenses that follow eq. (|4]|, 
shown by the dotted line, and form subsample SL0.5; solid open 
circles are for subsample SL0.75 (i.e. galaxies that follow eq. (5] 
shown by the solid lino). 



where (r, ip) are polar coordinates on the sky and 2 is the 
direction of the line-of-sight. Formally, the integral in eq. ^ 
should be calculated over —00 ^ z ^ +00, but we limited it 
over — lOr-cff ^ z ^ -l-lOr-cff. For an isothermal sphere this 
cut-off results in a negligible underestimation of the integral 
(« 4%). Coma galaxy velocity dispersions (Teff are measured 
by (1) reconstructing the line-of-sight velocity distributions 
(LOSVDs) from the kinematic moments, (2) coadding all 
LOSVDs inside r^s, each weighted by its projected light, 
and (3) fitting a Gaussian to the resulting LOSVD. They are 
listed in Tab.[T] Because of their higher signal-to-noise, only 
the major-axis data have been considered. The upper inte- 
gration limit rEin(o"eff ) in eq. (|6]) refers to eq. Q for the com- 
parison with SLACS subsample SL0.5 and to eq. ([SJ for the 
comparison with subsample SL0.75. In Fig. [5] we only show 
those SLACS lenses that belong to either subsample SL0.5 
(left-hand panels) or subsample SL0.75 (right-hand panels), 
respectively. 

The top row is for the total cylindrical mass of our 
standard two-component models with dark matter halos 
(i.e. with the best-fit density p = Tt^dyn x v + PDM.dyn in 
eq. [6} . The good agreement with the lensing results is re- 
assuring (the average mass offset is 0.05 dex for subsample 
SL0.5 and 0.02 dex for subsample SL0.75). It implies that 
the two completely independent methods yield consistent re- 
sults. Moreover, the scatter in the dynamical masses is not 
larger than in the lensing masses. Consequently, strong devi- 
ations from axisymmetry a r e unlik ely in the Coma galaxies. 
As shown in I Thomas et al.l l|2007ah . strong triaxiality, if not 
accounted for in the models, can bias dynamical masses by 
a factor of up to two, depending on viewing angle. Assum- 
ing random viewing angles, strongly triaxial mass distribu- 
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Figure 5. The projected total (luminous+dark) mass Meih within a fiducial Einstein radius rEin- Coma galaxies are indicated 
by the large symbols. Top row: two-component models with dark matter halos (p = T, dyn x J^ + PDM,dyn in equation|6ll; bottom 
row: dynamical mass under the as sumption that mass follows light (p = Tt.sc x i^)- Small circles: total projected masses of 
SLACS galaxies l|AuEer et al.ll2009l '). In the left-hand panels we compare Coma galaxies with the SLACS subsample SL0.5 using 
fiducial Einstein radii calculated via eq. Q. In the right-hand panels we compare to subsample SL0.75 (using eq. [5ll. Masses are 
plotted against the average velocity dispersion (Toft inside r^ff. 



tions would therefore likely cause a significant scatter in the 
dynamical masses, which is however not o bserved (but see 
Ivan de Ven. Mandelbaum. fc Keetonll2009l '). 

The bottom row of Fig. [5] is for projected masses of 
self-consistent dynamical models in which all the mass is as- 
sumed to follow the light (i.e. p = T,,sc x z/). These models 
are not consistent with the lensing masses. The discrepancy 
is larger for the SLACS subsample SL0.5 than for subsample 
SLo.75 because the smaller Einstein radii of subsample SL0.5 
emphasise the central r egions in the comparison. As it has 
been s hown earlier (e.g. iGerhard et al.|[200ll . iThomas et al.l 
l2007bl ). the mass distribution in early-type galaxies follows 
the light in the inner regions (pin ~ T,,dyn x Vin), but 
has an additional component in the outer parts (pout ~ 
T*,dyn X Vout + Pdm) ■ Then, assuming that all the mass fol- 
lows the light requires T,,sc > T,,dyn to include the outer 



dark matter (T*,sc x i^out ~ Pout ~ T,,dyn x i^out + Pdm). 
However, the central regions become proportionally more 

massive, too, such that Tt,ac X Uin > T,,dyii X Vin ~ Pin- 

This explains why the offset in the lower-left panel of Fig. [S] 
is larger than in the lower-right one. 



3.3 Stellar population masses 

In Fig. [6] we show projected stellar masses (i.e. using p = 
Tsaip X u and p — TKrou x v, respectively, in eq.[6|. The open 
dots in these panels represent the projected stellar masses 
for either the Salpeter IMF (top row) or the Kroupa IMF 
(bottom row). Kroupa stellar masses are always below lens- 
ing masses. The mass difference between the lenses and the 
Kroupa masses increases from subsample SL0.5 to subsample 
SLo.75 for the same reason discussed at the end of Sec. 13.21 
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Figure 6. As Fig. [5] but for stellar population masses. Top row: Salpeter IMF (p 

(p = TKrou X 1^)- 



Tsalp X u); bottom row: Kroupa IMF 



It could be due to dark matter. Salpeter stellar masses are 
likewise consistent with the lensing results. In low-dispersion 
galaxies the IMF cannot be much steeper than Salpeter 
as otherwise the implied stellar masses would exceed the 
total observed lens masses. In high-dispersion galaxies the 
Salpeter stellar masses are however not enough to explain 
the total lensing masses. Then, if all the lensing mass was 
stellar, the IMF would have to change with galaxy velocity 
dispersion and in massive early-types the stellar mass per 
stellar light would have to be larger than for a Salpeter IMF 
(or any equivalent top-heavy IMF). If the IMF is constant, 
then the top- left panel of Fig. [6] provides direct evidence for 
the presence of dark matter in high-dispersion early-type 
galaxies. 



3.4 Luminous and dark matter separated 

In our standard two-component models that take into ac- 
count the detailed, radially resolved stellar kinematics, the 
large projected masses of high-dispersion galaxies do not en- 



tirely originate from luminous mass. This is shown in the top 
panel of Fig. [7] (projected dark matter fractions inside the 
Ei nstein radius). N o te tha t here we compare to the results 
of iKoopmans et al.l (|2006l l. who provided a combined lens- 
ing and dynamics analysis of the first SLACS galaxies. We 
applied a similar lens selection as described in Sec. 13.11 (us- 
ing rEin = 0.03244 x acft —4.6324). Our dynamically derived 
dark matter fract i ons ar e in good agreement with those from 
iKoopmans et all ()2006h . As already stated above, they in- 
crease with (Toff. For comparison, we have also plotted the 
corresponding deprojected dark matter fractions (cf. eq. [3| 
inside the same three dimensional radius rEin(crefr) in the 
bottom panel of Fig. [71 They are generally lower and do 
not vary with a^g. Projection effects therefore contribute to 
the trends in Figs.[5][7l To get a better understanding about 
the stellar IMF and dark matter distribution it is necessary 
to analyse the intrinsic three dimensional properties of the 
galaxies. This will be done in the following Sees. 3] and [S] 
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4.1 Variation in the stellar IMF? 
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Figure 7. Dark matter fractions inside the Einstein radius. Top: 
projected dark matter fraction; bottom; deprojected dark matter 
fraction. Red/open circles: Coma galaxies of this work; dark/filled 
circles in the top panel: SLACS galaxies from 'Koop mans et akl 
||2006| '). Coma galaxy dark matter fractions are plotted against 
the velocity dispersion cTgff inside r^g. SLACS galaxies are plotted 
against the measured velocity dispersion inside the SDSS aper- 
ture. 



4 DYNAMICAL LUMINOUS MASS VERSUS 
STELLAR POPULATION MASS 

Fig. [8] shows dynamical and stellar-population mass-to-light 
ratios versus CTcff- The dynamical T»,dyn exhibit a broad 
distribution ranging from T,,dyn ~ 2 to T.^dyn ~ 10. In 
contrast, for the majority of Coma galaxies the stellar pop- 
ulation Tssp are almost constant with Tsaip « 5 — 6 for the 
Salpeter IMF and Txrou « 3 - 4 for the Kroupa IMF. In 
addition, while the dynamical T«,dyn clearly increase with 
galaxy velocity dispersion, there is no similar correlation be- 
tween (Tefi and stellar-population Tssp. 

There are two galaxies (GMP0756 and GMP1176) with 
distinctly lower stellar mass-to- light ratios than in the rest of 
the sample. As Fig. [9] shows, these two galaxies have overall 
young stellar populations. In general, one can read off from 
Figs. [5] and [5] that dynamical T^.dyn depend on a^g but not 
primarily on the stellar population age, while Tssp mostly 
reflect stellar ages and do not show any dependency on (Toff. 



Since dynamical and stellar population masses scale differ- 
ently with galaxy velocity dispersion (cf. Fig. the ratio 
Tt^dyn/Tssp has to vary with galaxy (Toflf (for any fixed IMF). 
This is explicitly shown by the large/red symbols in Fig. 1101 
In addition to our results, the figure also combines work from 
other groups. 

Pentagons represent the stellar- population analy s is of a 
subsample of our Coma galaxies bv lGrillo fc GobatI (|201(]| ) 
(cf. Sec. l2.2p . The pentagons differ from the large filled sym- 
bols only in terms of TKrou- The velocity dispersions and 
dynamical T,_dyn are the same. 

T riangles are for the SAURON survey l|Cappellari et all 
l2006h . In terms of both, the stellar population analysis 
(based on spectral absorption line indices) as well as the dy- 
namical modelling (orbit-based), they can be most directly 
com pared to the Com a galax ies of this work. Note, however, 
that ICappellari et al.] (|2006l ) measured only the sum of lu- 
minous and dark mass, assuming the latter to contribute 
only a small amount of mass in the central galaxy regions 
observed with SAURON (robs, max ^ ^eff). Adopting equiv- 
alent modelling assumptions for the Coma galaxies yields 
mass-to-light ratios typically 10 — 20 % higher than com- 
pared with stellar T*,dyn from models where dark matter is 
accounted for explicitly (cf. Tab. [T]). The overall distribu- 
tions of T,,dyii/TKrou are nevertheless similar in both sam- 
ples. 

Finally, Fig. [TU] also includes SLACS galaxies, analysed 
with a combined dynamics and lensing approach (Tr eu et al.l 
For these galaxies, the ratio of the total versus the 
stellar mass (the latter derived from broad-band colours) 
inside the Einstein radius (typically of the order of reff/2) is 
plotted along the y-axis. 

In all the samples included in Fig. [in] dynamical (or 
lensing) stellar masses systematically exceed the required 
masses for a Kroupa IMF. Above (Toff > 150 km/s the ratio 
T*,dyn/TKrou tcuds to iucrease with velocity dispersion. For 
lower mass galaxies the ratio becomes uncertain due to the 
more frequent presenc e of multiple stellar populations (e.g. 
ICappellari et al.ll2006l ). Moreover, for low-mass galaxies the 
assumption that all the mass follows the light is in conflict 
with gravitational lensing masses (cf. Sec.|3]). 

If the increase of T*,dyn/TKrou with a^a was a pure 
stellar population effect then we would have to assume that 
the IMF is not universal. Interpreted in this way. Fig. [10] 
would imply the IMF in high-dispersion galaxies to produce 
either more low-mass stars than the Kroupa IMF (i.e. be- 
ing Salpeter-like) or more stellar remnants (i.e. being top- 
heavy). A higher fraction of low- mass stars would reduce the 
number of SNe of type II (per stellar mass) and would lead to 
an overall lower metallicity. Thus, if the IMF changed from 
Kroupa towards Salpeter, then the increase in T«,dyn/TKrou 
would be expected to come along with a decrease in [Z/H]. 
Instead, a top-heavy IMF enhances the importance of type 
II SNe over type la SNe. Accordingly, if the IMF changed 
from Kroupa towards being top-heavy then one would ex- 
pect higher [g/Fe] in galaxies with higher T«.dvn/TKr ou 
^Thomas. Greggio. fc Bendeilll999l : Ig raves fc Fabedl2O10h . 

Fig. 111! shows Tt^dyn/TKrou agaiust stellar population 
age T, metallicity [Z/H] and [a/Fe] ratio. There is no corre- 
lation with any stellar population parameter. Note, however. 
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Figure 9. As Fig.|8l but stellar mass-to-light ratios are plotted against stellar population age 



that stellar metallicities and [a/Fe] ratios do not only de- 
pend on the stellar IMF, but also on the duration of the 
star- formation episode(s), the depth of the galaxy potential 
well and on evolutionary processes related to the cluster en- 
vironment. In this respect, the lack of evidence for an IMF 
change cannot be taken as a proof for a constant IMF. Al- 
ternatives to IMF variation are discussed in Sees. 14.41 and 

El 

4.2 The shape of the stellar IMF 

A direct comparison between dynamical and stellar- 
population mass-to-light ratios is provided by Fig. 1121 As 
already clear from the different scalings of T«,dyn on the 
one side and Tasp on the other (cf. Figs. [8] and [9|, neither 
the Kroupa IMF nor the Salpeter IMF yields a close match 
between T*,dyn and T 

ssp . 

Above the dotted line in Fig.[l2]the luminous dynamical 
mass is larger than the stellar mass required for the Kroupa 
IMF, while below the line the dynamical T,_dyn is formally 
insufficient for the Kroupa IMF. For the Salpeter IMF the 



corresponding limit is shifted towards T,,dyn which are a fac- 
tor 1.6 higher (dashed line). Accordingly, all Coma galaxies 
are compatible with a Kroupa IMF. The majority of the 
galaxies is also consistent with a Salpeter IMF, but there is 
at least one galaxy for which the dynamical T«,dyn is signif- 
icantly lower than Tsaip (at about the 3cr level; GMP5975). 
Concerning the total sample, however, the Salpeter IMF fits 
the dynamical masses better than the Kroupa IMF. The 
corresponding sample averages are (Tt^dyn/Tsaip) ~ 1-15 
for the Salpeter IMF (with an rms scatter of 0.35) and 
(T,,dyn/TKrou) ~ 1-8 for the Kroupa IMF, respectively. 

These conclusions also hold for other galaxy samples. 
Since Tsaip ~ 1.6 x Tkiou, the one-to-one line for the 
Salpeter IMF in Fig. 1101 would occur at T«_dyn/Taap ~ 1.6. 
Then, the Salpeter IMF provides on average a better match 
with dynamical/lensing masses. In line with this, recent 
near-infrared spectroscopic observations point towards a 
bottom-heavy IMF in massive earl y-type galaxies as well 
fvan Do kkum fc Conrovl[2010l . l201l h. Howeve r, our dynami- 
cal models as well as p revious lensing studies (jPerreras et al.l 
l2010l ; lTreu et al.ll201(l ) indicate that around tJeff < 200 km/s 
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models accounting for dark matter (this work); triangles: spectro- 
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matter jCappe llari et all 2006); pentagons: photometric Txrou 
l|Grillo fc Gobad l201(t and orbit-based dynamical models ac- 
counting for dark matter; open square s: photometric T xrou and 
combined lensing -|- dynamics models llTreu et al.ll201ol) . For all 
but the SLACS galaxies the average velocity dispersion (Teff inside 
red is plotted. The SLACS dispersions are the average over the 
spectroscopic aperture of the SDSS survey. 



the Salpeter stellar masses exceed the observed dynamical 
and/or lensing limits. This rules out a Salpeter IMF for low- 
mass galaxies. 

The results for the Coma galaxies are largely indepen- 
dent of the parameterisation chosen for the dark matter ha- 



los. Fits with logarithmic halos alone yield (T^'^^^^/Tsaip) = 
1.16, while NFW halos resuh in (T^_5^,/Tsaip) = 1.06. Both 
are consistent within the rms scatter (about ~ 0.32). 



4.3 Uncertainties in population Tasp 

Gas emission can refill the H/3 line and lead to an overesti- 
mate of stellar population ages and, then, of Tgsp. A young 
stellar subpopulation (dominating in terms of hght, but not 
in mass) can likewise bias stellar ages and Tasp, yet towards 
too low values. In any case, T,_dyn/Tssp would systemat- 
ically decrease with stellar population age. Fig. Illb shows 
however, that this is not the case in the Coma galaxies, such 
that a strong bias due to gas emission or young stellar sub- 
populations is unlikely. 

The stellar population parameters of the Coma galaxies 
are derived from spectral indices and, thus, represent aver- 
ages along the line-of-sight. In other words, at a given radius 
of observation robs, the SSP parameters combine the prop- 
erties of stars with r > r-obs- In contrast, dynamical mod- 
els are most sensitive to the mass distribution inside robs. 
Projection effects can therefore introduce a systematic bias 
between dynamical T,,dyn and stellar population Tssp if the 
stellar-population changes with radius. A monotonic stel- 
lar population gradient is enhanced after projection along 
the line-of-sight but diminished in the cum ulative mass-to - 
light ratio constrained by dynamical models (|Thomasll2006l ). 
More specifically, a radial increase of T leads to an underes- 
timation of T»,dyn/Tssp, a radial decrease of T to an over- 
estimation of Tt,dyn/Tssp. For a rather steep gradient of 
dlogT/dlogr = ±0.23 (a change by a factor of 1.7 per 
decade in radius), the expected systematic difference be- 
twee n T<..dvn and T ssp would amount to =f30 percent inside 
Vcs (| Thomas! |2006| ) . The observed gradients in the Coma 
galaxies are however much smaller (cf. Fig. [l|. Therefore 
systematics due to projection effects seem negligible. 
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Figure 12. Dynamical T, ^lyn versus stellar-population mass-to- 
light ratio Txrou- The dotted line shows the one-to-one relation 
for the Kroupa IMF. Since Ts^ip = 1.56 X Txrom the corre- 
sponding one-to-one relation for the Salpeter IMF occurs at larger 
(jyn (dashed line). 



4.4 Uncertainties and interpretation of dynamical 

The dynamical mass-to-light ratios Tt^dyn could be affected 
by systematic biases in the modelling process, e.g. arising 
from false symmetry assumptions. As it has been stated in 
Sec. [T] the luminous Tt.dyn can be biased by a factor of up 
to two, if the studied galaxies deviat e significantly from th e 
symmetry assumed in our models (| Thomas et al.l l2007al l . 
However, there is no evidence for the Coma galaxies to be 
strongly non-axisymmetric. Firstly, they do not show signif- 
icant isophotal twists as would be indicative for triaxiality. 
Secondly, as already discussed in Sec. [31 the good match be- 
tween dynamical and strong-lensing masses provides further 
evidence that the obtained dynamical masses are unbiased. 

Even if the luminous T*,dyn are accurate, they might 
not represent the galaxy stellar mass in a one-to-one fashion. 
Ambiguities can come from any non-stellar mass that follows 
the hght and contributes to T».dyn. 

Such a mass component could be gas loss during stellar 
evolution. It is not included in our Tssp, which only encom- 
pass the baryonic mass locked in stars or stellar remnants. 
For a 10 Gyr old population the lost gas mass amounts to 
about 40 perce nt of the originally formed stellar mass (e.g. 
lMarastonll2005h . Provided that it remains in the galaxies and 
provided it follows the hght distribution, it would contribute 
to the dynamical mass T,.dyn and the actual baryonic mass 
in stars and stellar remnants would only be (/*Tdyn) x L 
(with /, fa 0.6). Correcting the dynamical T*,dyn for stellar 
mass loss yields formally a good agreement with the Kroupa 
IMF: {(/*Tdyn)/TKrou> = 1.04 ±0.32. 

Theoretical arguments indicate that most of the lost 
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Figure 13. From top to bottom: total (luminous + dark) three 
dimensional mass density p, dark matter density PDM.dyn ^^nd 
luminous mass density p«_dyn- AH densities are spherically aver- 
aged. Left-hand panels: galaxies with less massive halos (ri^ ^ 
400 km/s); right-hand panels: galaxies with massive halos (djj > 
400 km/s. In the middle row the dotted lines indicate the spatial 
region where pDM.dyn < P*,dyn/10. 



gas is either expelled from the galaxies or recycled into new 
stars. Only a few percent of the original stellar mass is ex- 
pected to remain in hot gas halos a round the galaxies (e.g. 
ICiotti et aLlll99ll : lDavid et aLlliggi V Fittingly, observed gas 
mass fractions of hot X-ray halos around massive early-type 
galaxies are t ypically less than a percent of the present stel- 
lar mass (e.g. lMatsushitall200lD . such that stellar mass loss 
is unlikely to explain the excess T«,dyn > TKrou. 

A possible component of non-baryonic matter that fol- 
lows the hght is discussed below in Sec. 15.21 



5 DARK MATTER 

5.1 Mass that does not follow the light 

Fig. [13] shows the spherically averaged three dimensional 
density distributions of luminous and dark matter, as well 
as the sum of both. The sample is subdivided into galaxies 
with LOG-halo circular velocities Vh ^ 400 km/s and galax- 
ies with Vh > 400 km/s. The reason is that the latter galaxies 
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Figure 14. Comparison of the local logarithmic density slope ^ = d lnp(r)/d ln(r) with the density ratio PDM('")/p*,dyn(^) of 
dark to luminous matter at the same radius. Left panel: logarithmic slope of the luminous mass density p«,dyn = """t.dyn x 
right panel: logarithmic slope of the total mass density p = T, jyn x f + Pdm- Vertical dotted lines indicate where luminous 
and dark matter equalise. To the left of these lines luminous mass dominates (inner regions of the galaxies), to the right dark 
matter dominates (outer regions of galaxies). Horizontal lines are for f = —2.8, the luminous slope at which dark matter takes 
over luminous matter. 



have very uniform outer dark and total mass density profiles. 
The corresponding dark matter fractions (inside rcfi) scatter 
around a mean of {/dm, dyn) ~ 23± 17% and do not depend 
on (Tcff. Note that in the central galaxy regions our models 
implicitly maximise the mass contribution from the light. 
Similar maximum-bulge models for lensing galax ies yield 
dark matter fractions around 25 % as well (jBarnabe et al.l 
l2009l ). 



Fig. [14] shows the logarithmic slope of the luminosity 
density and the total mass density in the Coma galaxies. The 
slope is plotted against the ratio of dark to luminous matter 
densities at the same radii. In the inner regions the slope of 
the total mass density follows the light, since dark matter 
is negligible (p., dyn ^ Pdm), in the outer regions the total 
mass density profile is flatter than the light profile. Overall, 
the total mass density is roughly isothermal: p ~ r~^. This 
reflects t he nearly fiat circular velocity curves of ear ly-type 
galaxies jCerhard et all I2001I : pThomas et al.1 l2007bl ). Sim- 
ilar slopes for the t otal mass distribution have been seen 
in le nsing galaxies (|Koopmans et all I2OO6I : iBarnabe et al.l 
I2OIII ). Fig. [14] illustrates that the actual need for a dark 
matter component in our models comes from the fact that 
the outer mass distribution does not follow the light in early- 
type galaxies. The radius where the density of dark matter 
takes over luminous matter is roughly where the slope of 
the luminosity density falls below ^ ~ —2.8 (indicated by 
the horizontal dotted lines in Fig. I14p . 



5.2 A component of dark matter that follows the 
light? 

A galaxy might have more dark matter than captured by 
Pdm if some fraction of the halo mass follows the light so 
closely that it is mapped onto T,,dyn rather than pdm- In 
particular, Fig. 1 141 leaves the possibility open that this could 
happen in the inner galaxy regions where the slope of the 
luminosity distribution is —2 < ^ < —1. If the fraction of 
dark matter that follows the light is larger in galaxies with 
higher aea, then this would be a possible explanation for the 
trend between T,_dyn/TKrou and seen in Fig. [TU] 

Because the total (luminous + dark) mass Mtot,dyn is 
well constrained by the dynamical models (inside the region 
with kinematical data) , a spurious increase in the luminous 
mass component would be accompanied by a corresponding 
decrease in the nominal dark matter fraction /dm, dyn of the 
models. More specifically, under the assumption Mtot.dyn = 
Mtot,gai, the model parameters (T,,dyn, ./dm, dyn) and the ac- 
tual galaxy parameters (T,,gai, /DM,gai) would be related via 



Trt,dyn — T»,gal -|- (/DM,gal — /DM.dyn) Ttot.dyn, 



(7) 



where Ttot.dyn = Aftot.dyn/i is the total mass-to-light ratio 
(including dark matter). A degeneracy in the mass decom- 
position (at fixed total mass) would therefore correlate the 
offset 



AT = 



T,,dyn — T, 



gal 



(8) 



TKro u 

in stellar mass-to-light ratios with the offset 

A/dm = /DM.dyn " /DM,gal (9) 

in dark matter fractions as 
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Figure 15. (jy^/TKrou against dark matter fraction /DM.dyn (panel a) and against halo core radius rj, (scaled by the effective 
radius r^g ; panel b) . Shifting mass from the dark halo component into the luminous one at constant total mass moves galaxies 
along the direction indicated by the arrow in the left-hand panel (a). Note that for the right-hand panel the galaxy GMP1990 
has been omitted, since its dark matter fraction is so low that the determination of a halo core-radius becomes meaningless. 



AT : 



T 



tot.dyn 



T 



Krou 



A/d 



(10) 



Fig. 115b shows Tt.dyn/TKrou against the models' dark 
matter fractions /DM,dyn (inside rofi). There is a slight trend 
for T,_dyn/TKrou to be particularly large whenever /DM,dyn 
is low. Note that an intrinsic dark matter variation from 
galaxy to galaxy at a constant IMF scatters galaxies hori- 
zontally in Fig. 115b . while an IMF variation at constant dark 
matter fraction scatters galaxies vertically. The luminous- 
dark matter degeneracy discussed above scatters galaxies 
along the arrow shown in Fig. 115b . It marks the direc- 
tion along which the fitted (T,,dyn, /DM,dyn) are expected 
to separate from the galaxies' (T,.gai, /oM.gai) according to 
equation \iQ\ and - for each galaxy - depends on the ra- 
tio Ttot.dyn/TKrou. For the arrow in Fig. 115b we have used 
the average (Ttot.dyn/TKrou) = 2.39 over the Coma sample. 
The distribution of the majority of Coma galaxies roughly 
follows the arrow, in particular below /DM.dyn ^ 0.3. 

Fig. 115b suggests, though does not unambiguously 
prove, that the scatter in T*,dyn/TKrou could reflect a de- 
generacy in the dynamical mass decomposition. A contami- 
nation of some T*_dyn with dark matter would also explain 
the trend seen in Fig. 115b . where T,,dyn/TKrou is plotted 
against the halo core-radius rn (in units of ros). The more 
the inner dark matter goes into T,_dyn, the less the model 
component pdm traces the actual inner dark matter of the 
galaxies. Instead, it only represents the outer parts of the 
dark matter halos. Correspondingly, one would expect rela- 
tively larger halo core-radii (with respect to rcff) whenever 
Tt. dvn is large compa r ed to Tssp - as seen in Fig. [TSb. 

iNapolitano et ahl (|201G| ) analysed a large set of early- 
type galaxies with two-component spherical Jeans models 
and find dynamical T*,dyn close to the Salpeter IMF when 
using collisionless halos from cosmological simulations with- 



out baryon contraction. However, their T«,dyn get closer to 
the Kroupa IMF when baryonic contraction is taken into ac- 
count. Baryonic contraction might therefore be one way to 
make luminous and dark matter distributions similar enough 
to explain the difference between T,,dyn and Txrou (see 
Sec. El). 



In our modelling approach, the halo parameters are al- 
lowed to vary freely, without being connected to results from 
cosmological simulations. On the one hand this ensures that 
baryonic contraction is implicitly taken into account: the 
best-fit models for more contracted halos are simply ex- 
pected to occur in a different region of parameter space. 
On the other hand, the steepest halo density profiles that 
we probed are those from cosmological simulations without 
baryon contraction (NFW halos; cf. Sec. I2.1|l . If the actual 
galaxy halo profiles are steeper, then the best-fit dynamical 
model would still be obtained by shifting some fraction of 
the inner dark mass into T«_dyn- 

More similar distributions of luminous and dark matter 
in some galaxies than in others could also reflect differences 
in their evolut ionary histories. For example, the co smological 
simulations of iNaab. Johansson, fc OstrikeJ (|2009 ) indicate 
a difference in the radial distribution of in-situ formed stars 
relative to stars that were accreted during mergers. In-situ 
formed stars have a more centrally concentrated radial dis- 
tribution than stars that were accreted in collisionless merg- 
ers. The latter dominate the stellar mass density around roff . 
In any case, more detailed investigations of numerical simu- 
lations are required to conclude about a possible degeneracy 
between luminous and dark matter in the inner regions of 
galaxies. 
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Figure 16. Dark matter fractions (inside tlie effective radius reft ) 
against velocity dispersion a^ff. Open circles: fraction of dark mat- 
ter /DM.dyn that does not follow the light; filled circles: dark mat- 
ter fraction /dm Krou assuming i) a universal Kroupa IMF and 
ii) that the excess mass (T,, jy^^ — TKrou) x L is a component of 
dark matter that follows the light. 




Figure 17. Spherically averaged dark matter density Pdm Krou 



against radius. The excess mass (T, 



dyn " 



-Txrou) xi with respect 



to a Kroupa IMF stellar population has been added to the dark 
matter halo. The dashed line indicates a power-law density with 
logarithmic slope of -2. 



5.3 The distribution of dark matter in case of a 
universal Kroupa IMF 

If we adopt the point of view that the stellar IMF in early- 
type galaxies is universal and Kroupa-like then this affects 
the distribution of dark matter significantly. The reason is 
that in this case our nominal halo component captures only 
a part of the galaxies' dark matter, while a large fraction 
follows the light and is included in T,,dyn- In fact, a universal 
Kroupa IMF implies the dark matter fractions /dm, Krou to 
read 

J- / \ — s / \ , ('""•.dyn ~ TKrou) X L[r) 

/DM,Krou(r-j = /DM.dyn (r) H — -— . (11) 

AAot,dyn(r) 

These fractions are larger than the nominal /dm, dyn derived 
from Pdm, have a smaller scatter and slightly increase with 
galaxy a^s (cf Fig. [T6)) . 

Fig- El shows the spherically averaged dark matter den- 
sity profiles 



pDM.Krou = pDM + (T,_dyn — TKrou) X U 



(12) 



including the extra dark matter required for a Kroupa IMF 
(cf. equation [2}. These density profiles are smooth and close 
to a power-law with logarithmic slope slightly shallower than 
-2. The slight wiggles in the profiles around ~ 5 kpc indicate 
the transition from the outer parts, which are dominated by 
Pdm, to the inner parts, which are dominated by the second 
term on the right-hand side of equation (|12l) . 



5.4 Dark matter density and halo assembly epoch 

In lThomas et al.l (l2009bh we estimated halo assembly epochs 
Zform based on the assumption that the average dark mat- 



ter density (pdm) inside 2 reft scales with (l-l-ztorm) • From 
the overdensity of dark matter in early-type relative to spi- 
ral galaxies one can then narrow down elliptical galaxy as- 
sembly redshifts with an additional assumption about the 
typical formation redshift of spiral galaxies (zform ~ !)• 
Adopting a universal Kroupa IMF results in average dark 
matter densities a factor of ~ 3 larger than the nominal 
ones of the dynamical models. The strongest contribution to 
this increase comes from the single galaxy GMP1990 which 
has a negligible /DM.dyn, but a large T*,dyn/TKrou. Without 
this galaxy, dark matter densities are only larger by about 
a factor of 1.6, which lets the halo assembly re dshifts in- 
crease from Zform ~ 1 — 3 (|Thomas et al. I l2009bl ) to about 

Zforia ^ 1-5 3.5. 

As a second test, we have compared the poM.Krou 
from Fig. 1171 direct l y to the galaxy formation models of 
Lucia fc Blaizod i|2007l ). These models do not include 
the dynamical reaction of dark matter on the baryon infall. 
Therefore, we first subtracted from the observed dark mat- 
ter profile pDM.Krou the expected effect of baryon contraction 
by assuming the adiabatic approximation (i.e. by inverting 
the equations of Blumonthal ct al. 198^). The relationship 
between the average dark m atter density and ha l o form a- 
tion redshift in the models of iDe Lucia fc Blaizod (|2007l ) is 
well fitted by log(pDM) ~ —2.9 -I- log(l -|- Zform)^- We used 
this relation to calculate Coma galaxy assembly redshifts 
from the decontracted Kroupa dark matter halo densities. 
Fig. llSl shows that these halo assembly redshifts cover a sim- 
ilar range (zform « 1 — 3) as the ones obtained through the 
comparison with spiral galaxy dark matter densities. Note 
that we here computed star-formation re dshifts from the av - 



erage stellar ages inside TcB, while in iThomas et al 



we used the central stellar ages from lMehlert et al 
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Figure 18. Star formation epoch (x-axis) against halo assembly 
epoch (y-axis). The dotted line represents the one-to-one relation. 



For many of the Coma galaxies, star formation red- 
shifts and dark halo assembly redshifts are similar. Two 
galaxies appear far to the left of the one-to-one relation. 
These are GMP0756 and GMP1176, which have extended 
and relatively young (r ~ 3 Gyr) stellar disks. The galaxy 
on the bottom-right, where the stars seem sig nificantly older 
than the halo is GMP5568. As discussed in omas et al.l 
(|2009bl ). dry mergers can move galaxies below the one-to- 
one relation in Fig. [181 Besides this possibility the effective 
radius of GMP5568 is exceptionally large [r^s ~ 27kpc). 
Therefore, the dark matter density refers to a spatial region 
significantly larger than in any of our other Coma galax- 
ies. If the effective radius of this galaxy was overestimated, 
then the average dark matter density would be underesti- 
mated and so would be the halo assembly redshift. This 
could also explain the high dark matter fraction of GMP5568 
(cf Tab.[l|, since it generally increases with the physical dis- 
tance from the galaxy centre. 



6 THE TILT OF THE FUNDAMENTAL PLANE 

The effective radius , mean surface brightness (/)cfi inside 
rcff and the central velocity dispersion ao of a galaxy are 
connected via 

CM (13) 



2 



'■off 



where M is the total mass and 
L 



2tt r\ 

cri 



(14) 



where L is the total luminosity. In virial equilibrium the 
structure coefficient cm depends on the orbital structure and 
the radial distributions of the mass and the tracer popula- 
tion. For a homologous family of dynamical objects cm is a 



constant and with (T) = M/L equations (|13p and (|14|) lead 
to 



lo ■ — = 2 lo 

kpc km/s 

and 



log T ^,,9 + 7 



Lqpc- 



log 



(T> 



M0/Lfl,c 



log 2tt 



(km/s)2kpcM5 



rCM 



(15) 



(16) 



The actually observed fundamental plane (FP; 
iDiorgovs ki fc Davisiri987l : [Dressier et al.lll987h of early-type 
galaxies reads 



log 



kpc 



Q log 



km/s 



+ P log 



We 



^QPC- 



+ 7 



(17) 



with a ^ 2 and /3 7^ —1. It is tilted with respect to the virial 
plane of equation psp . 

Fig. 119b shows the FF0 of Coma galaxies. The best-fit 
parameters of an orthogonal fit are a = 1.26 ± 0.58 and 
P = —0.72 ± 0.10 (bootstrap errors). Within the statisti- 
cal uncertainties, the fit is consistent with the r-band FP of 
SPSS early-type galaxies (a = 1.49±0 05, 13 = -0.75±0.01; 
iBernardi et~aLll2003l : iHvde fc Bernardill2009l ). The larger er- 
rors result from the smaller sample size. Most of the differ- 
ence with respect to Bcrnardi ct al. ( 2003i) goes back to four 
galaxies (GMP0144, GMP0756, GMP1176, and GMP5975) 
that are distinct from the rest of the sample in several re- 
sp ects: (1) they have y oung central stellar cores (ro < 7 Gyr; 
cf. iMehlert et al]|2003l ): (2) at least two of them have an ex- 
tended thin stell ar disk; (3) they follo w different dark halo 
scaling relations l|Thomas et al.ll2009a) . Fig.l^Hlis equivalent 
to Fig. ll9l except that these four galaxies have been removed. 
The corresponding FP matches well with the SDSS results. 



6.1 The fundamental mass plane 

The tilt in the FP can reflect non-homology (i.e. cm being 
a function of galaxy mass) or variations in (T) (or both) . A 
systematic variation of (T), in turn, could reflect stellar pop- 
ulation effects or changes in the dark matter distribution. In 
any case, the dynamical models allow to incorporate varia- 
tions of (T) into the FP. For this purpose let (S)cff denote 
the average surface mass density inside Tad. Then, equation 
P7|l can be written 



log 



Tcff 

kpc 
where 

7 = - lO! 



a log 



km/s 



+ /3 1og 



M, 



0pc- 



+ 7, 



G 



(km/s) 2 kpc 



-Cm 



(18) 



(19) 



Equation ([IS t defin es the so-called fundamental mass plane 
iBolton et al.l l2007'). 

Figs. 119b - I19H show Coma galaxy MPs for different 
choices of {S)cfi . Using either (E)Krou = Tktou x (/)eff (panel 
b) or {E),_dyn = T,_dyn X (/)cff (panel c), the tilt is re- 
duced, but does not vanish (/3 ~1)- The change in a 



2 Note that in this section we use the average velocity dispersion 
(70 inside the central 2" , derived in the same way as the effective 



cTcff discussed in Sec. \3\ 
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Figure 19. Fundamental plane and mass plane of 16 Coma galaxies. In each panel the best-fit parameters of an orthogonal 
fit and the rms-scatter in log r^ff are quoted, a) Traditional FP; b) fundamental mass plane from stellar populations, i.e. the 
effective surface brightness is (/)efi is replaced by the mass-density (S)efi = TKrou x (^)eflFi c) as b) but for the dynamically 
derived stellar-population T,_dyii: {S),_dyn = '^*,dyn X (^)cff; d) as c), but instead of the dynamically derived steiiar mass-to- 
light ratio the total Mtot.dyn/-^' (including dark matter) is used: (S)tot_(jyn = A/tot, dyn/-^^ x (^)Gft (-'^^tot,dyn/-^^ is taken at the 
effective radius). Solid lines trace the one-to-one relation. 



from F ig. 119b. to Fig. | 19b i s con sistent with the SDSS anal- 
ysis of iHvde fc Bernardil (|2009l ). Fig. [T9h shows the case 
of (E)efr including all the projected mass (luminous and 
dark) inside rcff , respectively. The mass plane of Fig. I19t i 
is tilt-free within the errors. The match to the virial plane 
of equation (|15p improves further when omitting the four 
galaxies harbouring young stellar cores (cf. Fig. I20|) . With 
(T) = Mtot,dyn/i the MP has the same scatter as the FP 
itself (cf. Figs. 119b and d). For the subsample of old Coma 
galaxies the scatter slightly increases (cf. Fig. I20p . 



from TasjEI. Th e amount of reduction is consistent with the 
SPSS results oflHyde fc Bernardil (|2009l ). In accordance with 
[Graves fc Fabei (|2010l ). the tilt is further reduced if Tssp 
is replaced by the dynamical Tt,dyn- As discussed above, 
the different scalings of T,,dyn and Tssp with acfi could re- 
flect (1) changes in the IMF or (2) changes in the distribu- 
tion of dark matter. The absence of any correlation between 
T*.dyn/TKrou and stellar population parameters makes (2) 
more likely than (1). Finally, the tilt vanishes completely 
(for a subsample of Coma galaxies with uniformly old stellar 
populations), if the remaining dark mass inside the effective 



6.2 The tilt of the fundamental plane 

The tilt in the FP is reduced when the effective surface 
brightness is replaced by an effective surface mass derived 



In Sec. \E\ the surface mass was derived from the Kroupa IMF, 
but the difference between the Kroupa and the Salpcter IMF is 
only a constant scaling factor. The results for the Salpeter IMF 
are therefore similar. 
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Figure 20. As Fig. 1191 but four galaxies with young stellar cores are omitted. 



radius r^s is taken into acco unt as well. This b ehaviour is 
also found in lensing studies l|Bolton et al.|[2007l ). 

In our FP analysis we have not tried to calculate cm 
directly from the dynamical models. However, the density 
distribution (in particular the galaxy flattening) as well 
as the orbital struct ure vary among the Coma galaxies 
l|Thomas et al.ll2009al '). Nevertheless, Fig.[20ji indicates that 
for old Coma early-types the tilt of the FP is dominated by 
mass-to-light ratio effects rather than any possible variation 
in Cm- 

6.3 Mass estimators 

IWolf et all l|20ld ) provide an estimation for the mass 

Mwio(r3) = ^^^^ (20) 

inside the radius 7-3 where the logarithmic density slope of 
the tracer population is ^ = —3. Here, (cr^) i s the aver- 
age of the projected a over the whole galaxy. IWolf et al.l 
l|20ld ) derived eq. (O for non-rotating spherical galaxies. 



For the Coma galaxies, we averaged the measured up 
to Tcff , ignoring the galaxies' rotation velocities. Fig. 121b 
compares Coma galaxy masses Mtot,dyn(7'3) against the pre- 
dictions of eq. (|20p . Averaged over all Coma galaxies we 
find (Mtot,dyn(r3)/Mwio(r3)) = 1.03 ± 0.27 (rms scatter). 
Fig. 121b is similar to Fig. 121b ,. except for the additional 
approximation ~ ri/2 ~ 4/3 ros, where ri/2 is the de- 
projected half-light radius. The agreement with the Coma 
galaxies is still very good. 

Fig. [22] compares luminous dynamical masses with the 
estimator 

Mcoe = ^ . (21) 

from lCappellari et al.1 (|2006l ). When using luminous dynam- 
ical masses M,,dyn = T,,dyn X L from model fits that 
do have a separate dark matter component, then we find 
(M*,dyn/Afco6) = 0.86 ±0.35 (cf. Fig. [Ha,). The small offset 
disappears if we use model fits that do not have an addi- 
tional dark mat ter component (equiva lent to the approxi- 
mation made in ICappellari et aUl2006r ). The corresponding 
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Figur e 21. Comparison of Coma galaxies witli the IWoIf et alj 
mass estimator (cf . eq. I20I I : (a) total dynamical mass en- 
closed inside the radius ra, where the logarithmic slope of the 
luminosity distribution equals § = —3; (b) is similar to (a) but 
with the additional approximation ra fa i/^r^ff. 



Figure 22. C ompa rison of Coma galaxies with the 
ICappellari et aU 1120061) virial mass estimator (cf. eq. I21II : 
(a) luminous dynamical mass X L from fits explicitly 

allowing for a separate component of dark matter; (b) luminous 
dynamical mass T, sc X L from fits assuming that all the mass 
follows the light. 



luminous dynamical masses Mt,Bc = T*,sc x L (cf. Tab. [T]) 
are slightly larger and the comparison with eq. (I2ip yields 
(M.,sc/Mco6) = 1.01 ±0.36 (cf. Fig. [lU. The rms scatter 
include s both the measureme nt errors and anisotropy vari- 
ations (Tho mas et al ]|2007b"). Note, however, that the as- 
sumption that all the mass follows the hght is inconsistent 
with lensing m asses (cf. Fi g . [5]). 

Thus, the IWolf et all (|2010l l formula gives good esti- 
mates for the total dynamical mass inside a radius which is 
a bit l arger than r^a- The virial estimator of jCappellari et al] 
l|2006l ) captures the entire dynamical mass that results un- 
der the assumption that total mass follows light. Assuming 
rd. ~rii2, eq. ([21} implies 

u A/f I \ ^C06 2.5a-gffreff , . 

Mcoelrs) ~ Mco6(ri/2) = — ^ — = , (22) 

whereas from eq. 1201 and ~ 4/3roff it follows 

Mwio(r3) ^ 11^^. (23) 

Concerning the enclosed mass inside 4/3rcff, the two mass 
estimators would differ by a factor 4/2.5 = 1.6, unless 
{a^) 7^ cr^fj. Due to the different treatment of rotation 
CTgfj « 1.3 (cr^) in the Coma sample, such that the actual 
difference is only about 20 percent. 



7 SUMMARY 

We compared dynamically derived stellar mass-to-light ra- 
tios Tt^dyn with completely independent results from simple 
stellar population models. Our dynamical models are based 
on Schwarzschild's orbit superposition technique and have 



two mass components. One follows the light and its mass- 
to-light ratio T»,dyn is assumed to approximate the stellar 
mass distribution. The other mass component explicitly ac- 
counts for dark matter. This way, any potential degeneracy 
between the stellar mass and the dark matter halo is min- 
imised. The Coma galaxy sample studied here is currently 
the largest with axisymmetric Schwarzschild models includ- 
ing dark matter explicitly. 

Intrinsic uncertainties in the modelling, in particular re- 
lated to the assumption of axial symmetry, are unlikely to 
bias our results significantly. The main reason is that in pro- 
jection, our dynamical masses match well with completely 
independent results from strong gravitational lensing. 

Our main findings are: 

(i) For galaxies with low velocity dispersions (tTofi ~ 
200 km/s), the assumption that all the mass follows the light 
yields projected masses larger than in comparable strong- 
gravitational lens systems. 

(ii) In high- velocity dispersion galaxies (acff ~ 300 km/s) 
the assumption that mass follows light is consistent with 
strong lensing results. 

(iii) Two-component dynamical models with an explicit 
dark halo component yield total projected masses that are 
in good agreement with results from strong gravitational 
lensing for all galaxies. 

(iv) In two-component models, the mass-to-light ratio 
T,,dyn of the component that follows the light increases with 
galaxy velocity dispersion a^s- 

(v) Stellar population Tssp (for any fixed IMF) are largely 
independent of OcS- As a result, the ratio T,^dyn/Tssp of 
luminous dynamical mass over stellar population mass in- 
creases with galaxy velocity dispersion. 
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(vi) The luminous dynamical T,,dyn is always larger than, 
or at least equalises, the stellar-population mass-to-light ra- 
tio Tkfou for a Kroupa IMF. 

(vii) There is no correlation between Tt_dyn/"''"ssp ^'^'^ 
stellar population age, metallicity or [a/Fe] ratio. 

(viii) Inside r^s, the average fraction of dark matter (that 
does not follow the light) is /oM.dyn = 28± 17 % in the Coma 
galaxies. 

(ix) The tilt of the FP reduces if the effective surface 
brightness (/)Gff is replaced by the stellar population surface- 
mass density Tssp x {/)cff , further reduced if (/)cft is replaced 
by the dynamical stellar surface-mass density T,,dyn x {/)cff 
and, for a subsample of galaxies with uniformly old stellar 
populations, vanishes completely with the total dynamical 
surface-mass density (Stot,dyn)cfi. 

(x) Commonly used mass estimators are accurate to the 
20 - 30 % level. 

The implications of these findings are as follows; 

(i) That luminous dynamical masses increase more 
rapidly with galaxy velocity dispersion than stellar- 
population masses for a fixed IMF could be due to a change 
in the IMF or due to an increasing amount of dark matter 
following a spatial distribution similar to that of the light. 

(ii) If the IMF changes, then massive early-types (tTcff ~ 
300km/s) have up to two times more stellar mass per stellar 
light than lower-mass galaxies ((Joff ~ 200km/s), which are 
consistent with a Kroupa IMF. However, the lack of any 
correlation between T«,dyn/Taap and stellar population age, 
metallicity or [a/Fe] ratio is consistent with, though does 
not prove, that the IMF is actually universal. 

(iii) If the IMF is universal, then the increase in luminous 
dynamical masses must primarily come from a component 
of dark matter that follows the light very closely and is more 
important in more massive galaxies. The IMF would be con- 
sistent with being Kroupa in all early-types. 

(iv) Independent of the actual slope of the stellar IMF, 
luminous dynamical masses are on average more accurately 
predicted by assuming a Salpeter IMF: (T,^dyn/Tsaip) ~ 
1.15, but these masses may not represent exclusively stars. 
The Kroupa IMF yields (T.,dyn/TKrou> = 1.8. 

(v) Adopting a Kroupa IMF and counting the excess mass 
(Tt,dyn — TKrou) X L as dark matter that follows the light 
doubles the average dark matter fractions inside r-off to about 
/DM,Krou ~ 55±12 %. Moreover, it yields a smooth trend be- 
tween the resulting /DM,Krou and galaxy velocity dispersion 
and, also, smooth dark matter halo profiles. 

(vi) The FP tilt is not a pure stellar population effect. 
Further inferences about the tilt depend on the interpreta- 
tion of the observed T»,dyn/TKrou. As above, that the tilt 
reduces when considering the dynamical mass Tt_dyn x L 
that follows the light could be due to (1) variations in the 
relative distribution of luminous and dark matter or (2) IMF 
variability. 
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APPENDIX A: ERRATUM: DYNAMICAL 
MASSES OF EARLY- TYPE GALAXIES: A 
COMPARISON TO LENSING RESULTS AND 
IMPLICATIONS FOR THE STELLAR INITIAL 
MASS FUNCTION AND THE DISTRIBUTION 
OF DARK MATTER 

The paper 'Dynamical masses of early-type galaxies: a com- 
parison to lensing results and implications for the stellar 
initial mass function and the distribution of dark matter' 
was published in Mon. Not. R. Astron. Soc. 415, 545-562 
(2011). A corrected version of fig. 5 of that paper is given 
here along wit h affected statements in the text. 

In fig. 5 of'Thomas et all (|201ll ') we compared projected 
masses of Coma galaxies obtained (1) from dynamical mod- 
els with explicit dark matter halos (top row) and (2) from 
dynamical models in which all the mass follows the light 
(bot tom row) with grav itational lenses from the SLAGS sur- 
vey (|Auger et a l."2009^'). The Coma galaxies' mass estimates 
without dark matter (open circles in the bottom panels) were 
inadvertently raised by 0.19 dex. Fig. I All shows the correct 
values. The Coma galaxies' masses of our fiducial models 
with dark matter, shown in the upper panels, were correct 
in the original paper version and we repeat the upper panels 
for the sake of completeness. 

Neglecting outer dark matter in dynamical models 
where all the mass is assumed to follow the light leads to 
a mass deficit. Therefore, the masses obtained from these 
models fall below those of measured strong gravitational 
lenses. The effect becomes noticeable at larger radii where 
the luminous mass is less important. This can be seen in the 
bottom panels of Fig. lAll the left-hand panel is for lenses 
which have an Einstein radius rEin ~ 0.5 r^ff whereas the 
right-hand panel is for lenses with rBin ~ 0.75 reff. At 0.5 retr 
projected masses of dynamical models without dark mat- 
ter are still consistent with the lenses, while at 0.75 rcff the 
mass deficit in the outer parts becomes apparent. Dynami- 
cal models with dark matter (upper panels) are consistent 
with gravitational lens masses at all radii. 

Con sequently, the d escrip tion of the bottom panels in 
fig. 5 of iThomas et al] l|201ll ). last paragraph of Sec. 3.2, 
'The bottom row of Fig. 5 [...] than in the lower-right one.' is 
incorrect, while the conclusion drawn from the comparison - 
that dynamical models in which all the mass follows the light 
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are inconsistent with strong gravitational lensing results - 
remains correct. 

In addition, two statements in the Summary (Section 7) 
are incorrect: '(i) For galaxies with low velocity dispersions 
((Tcff ~ 200km/s), the assumption that all the mass follows 
the light yields projected masses larger than in comparable 
strong-gravitational lens systems.' and '(ii) In high- velocity 
dispersion galaxies (cTefi ~ 300km/s) the assumption that 
mass follows light is consistent with strong lensing results.'. 
Instead, the discrepancy shows up as a mass deficit and does 
not depend on the galaxy velocity dispersion CToff- 

We note that in the rest of the paper we only discuss 
the properties of our fiducial dynamical models with dark 
matter and the conclusions are unaffected by this mistake. 

This paper has been typeset from a Tj^X/ I^Ti^X file prepared 
by the author. 
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Figure Al. Corrected version of fig. 5 from iThomaa et al.l ll201lh . Top panels (as in the original paper): the projected total 
(luminous+dark) mass Meih within a fiducial Einstein radius rBin from two-component dynamical models with dark matter 
haloe s. Coma galaxies are indicated by the large symbols, while the small circles are SLACS gravitational lenses jAuger et al.l 
I2OO9I V Bottom panels: similar projected mass, but from dynamical models in which all the mass follows the light. I n the left-hand 
panels the comparison is made at rgin ~ 0.5 r^g, while in the right-hand panels at rEin ~ 0.75 r^g (details in lThomas et akl 
Eolll 'l. In the original version of the paper the Coma galaxies' masses in the bottom panels (open circles) were erroneously raised 
by 0.19 dex. 
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